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INTRODUCTION 


Nucleate  pool  boiling  is  a  very  efficient  heat 
transfer  process  for  removing  heat  from  a  heat  transfer 
surface.  It  is  one  of  the  important  mechanisms  to  remove 
heat  from  electronic  and  microelectronic  equipment, 
nuclear  reactors,  space  vehicles,  and  other  heat  transfer 
equipment  used  in  process,  refrigeration,  and  food 
industry.  In  order  to  understand  the  process  of  heat 
transfer  it  is  important  to  understand  the  underlying 
mechanism  of  the  process.  An  analytical  model  consistent 
with  the  requirements  of  nucleate  pool  boiling  heat 
transfer  requires  mathematical  expressions  for: 

1 )  Number  of  Nucleation  Sites 

2)  Bubble  Departure  Diameters 

3)  Bubble  Growth  Rates 

4)  Bubble  Emission  Frequency 

The  purpose  of  this  work  is  to  calculate  the  bubble 
departure  diameters  and  bubble  emission  frequency  for 
the  data  of  Engelhorn.  Engelhorn  has  conducted 
experiments  for  a  wide  variety  of  refrigerants  and  for  a 
wide  range  of  heat  flux  and  pressure.  For  the  calculation 
of  bubble  departure  diameters  we  used  Laplace  Equation 
and  for  the  bubble  emission  frequency  we  used  the 
equations  earlier  developed  by  Sharma  et  al. 
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ANALYSIS: 


Bubble  Departure  Diameter 


The  bubble  departure  diameters  are  obtained  by  analyzing 
the  forces  on  a  typical  growing  bubble.  Two  forces  are 
important;  surface  tension  and  the  buoyancy. 

Surface  Tension  Force  =  tc  D  a 

Buoyancy  Force  =  [pj  -  p^]  g 

At  equilibrium,  the  two  forces  are  equal. 

^  D3  [pj  -  py]  g  =  71 D  a 

^Vlpi  -"pvte 

Equation  (1)  is  known  as  the  Laplace  Equation.  The 
equation  shows  that  bubble  departure  diameters  is  the 
function  only  of  physical  properties  which  are  constant  for 
a  given  pressure. 
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Bubble  Emission  Frequency 

In  nucleate  pool  boiling  heat  transfer,  bubble  emission 
frequency  (f)  is  composed  of  two  periods; 

1)  Growth  period  (0d)  and  2)  Waiting  period  (0w) 

Therefore,  bubble  emission  frequency  f  is: 

f= - ^ 

(0d  +  0w) 


The  equations  were  earlier  developed  for  the  growth 
period  and  the  waiting  period.  From  those  equations,  the 
final  equations  for  the  bubble  emission  frequency  are: 


for  Ja  <  100 
f  = 


fl33.3/Pl2[g/(pi  -  pv)g  0.867 
It  al  Ja2  al 


qw 


and  for  Ja  >100 
f  = 


[133.3/P]2[g/(pi  -  pv)g 

25alJa3/2  + 


0.867  pkjAT w^2 
ocl  qw  ^ 


The  above  equations  reveal  that  for  a  given  liquid  and  at  a 
given  pressure  the  bubble  emission  frequency  is  a 
function  of  heat  flux,  and  physico-thermal  properties  . 
These  equations  were  used  to  calculate  the  bubble 
emission  frequency  for  the  data  of  Engelhorn. 
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RESULTS 


The  bubble  emission  frequency  were  calculated  for 
refrigerants  R-11,  R-12,  R-13,  R-13B1,  and  R-22. 
Table- 1  shows  the  values  of  pressures  for  each  fluid  and 
the  range  of  heat  flux. 


TABLE- 1 


REFRIGERANTS 

HEAT  FLUX, 

W/m2 

PRESSURE,  Bar 

R-11 

1,000  to  83,000 

0.019,0.028, 
0.503,  &  0.991 

R-12 

100  to  102,000 

0.25,  0.50, 

1.0  &1.80 

R-13 

200  to  84,000 

2.80,  4.55, 

7.35  &10.55 

R-13B1 

200  to  96,000 

0.78,1.51* 

5.60 

R-22 

200  to  99,000 

0.39,0.84  & 

2.15 

Graphs  were  prepared  for  the  frequency  as  a  function  of 


heat  flux.  The  study  reveals  that  bubble  frequency  is  the 
strong  function  of  heat  flux.  It  increases  with  increase  in 
heat  flux.  The  inspection  of  these  graphs  also  reveal  that 
bubble  frequency  also  is  a  function  of  pressure.  An 
increase  in  pressure  shows  increase  in  bubble  emission 
frequency.  An  increase  in  bubble  emission  frequency 
with  increase  in  heat  flux  and  pressure  implies  that  heat 
transfer  rates  in  boiling  will  be  enhanced  with  increase  in 
heat  flux  and  pressure. 


The  calculations  also  reflect  the  effect  of  pressure  on 
bubble  departure  diameter.  The  bubble  departure 
diameters  decrease  with  increase  in  pressure.  This 
functional  relationship  is  shown  by  making  plots  of  bubble 
departure  diameter  as  a  function  of  pressure  for  five 
refrigerants  investigated.  The  increase  in  the  bubble 
departure  diameter  with  decrease  in  pressure  will 
manifest  as  reduction  in  heat  transfer  coefficients  as 
pressure  decreases. 


CONCLUSIONS 


1.  The  study  reveals  that  the  bubble  emission 
frequency  is  a  strong  function  of  heat  flux.  It  also 
increases  with  increase  in  pressure. 

2.  Bubble  departure  diameters  are  the  function  only  of 
pressure  for  a  given  liquid.  The  bubble  departure 
diameter  decreases  with  increase  in  pressure. 
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Power  Curve  Fit 

Pts  Plotted  =  4  Offscale  Pts  =  0 

Regression  Equation: 

Y  =  9.80831  E-09  X  1.90664 

Corretation  Coefficient  *  .998636 

X-axis  file:  r11_heatflux_1 
Y-axis  file:  r11_freq_1 
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Power  Curve  Fit 
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Regression  Equation: 
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f  VS  q  at  0.503  bar  for  r-1 1  (Engelhorn) 
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Power  Curve  Fit 

Pts  Plotted  =  8  Offscale  Pts  =  0 

Regression  Equation: 

Y  =  1 .7335E-04  X  1 .08382 

Correlation  Coefficient  =  .990954 
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Power  Curve  Fit 

Pts  Plotted  =  11  Offscale  Pts  =  0 

Regression  Equation: 
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X-axis  file:  r12_heatflux_1 
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Logarithmic  Curve  Fit 

Pts  Plotted  =  8  Offscale  Pts  =  0 

Regression  Equation: 

Y  =-49.7848  +  (  6.38598  )  LNX 

Correlation  Coefficient  =  .984358 


X-axis  file:  r12_heatfkjx 
Y-axis  file:  r12_freq 
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Logarithmic  Curve  Fit 

Pts  Plotted  =  5  Offscaie  Pts  =  0 


Regression  Equation: 
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Correlation  Coefficient  =  .993208 

X-axis  file:  r12_heatflux2 
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f  vs  q  at  2.8  for  R-13  for  the  Data  of  Engelhorn 
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Power  Curve  Fit 

Pts  Plotted  =  9  Offscale  Pts  =  0 

Regression  Equation: 

Y  =  1.15874E-03  X  .93721 

Correlation  Coefficient  =  .99093 

X-axis  file:  r13_beatflux_1 
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f  VS  q  at  4.55  bar  for  R-13  (Engelhorn) 


Heat  Flux  w/sq.  m 

Report  Created:  04-17-1996  10:13:31  AM 

Logarithmic  Curve  Fit 

Pts  Plotted  =  5  Offscaie  Pts  =  0 

Regression  Equation: 

Y  =-31.0799  +  (4.20439)  LNX 

Correlation  Coefficient  =  .981057 

X-axis  file:  r12_beatflux1 
Y-axis  file:  r13_freq1 
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Regression  Equation: 
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f  VS  q  at  10.55  bar  for  R-13  (Engelhorn) 


Power  Curve  Fit 
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Regression  Equation: 
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Power  Curve  Fit 

Pts  Plotted  =  9  Offscale  Pts  =  0 

Regression  Equation: 

Y  =  3.8682E-04  X  .999344 

Correlation  Coefficient  *  .99372 

X-axis  file:  r22_heatflux_1 
Y-axis  file:  r22_freq_1 
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f  VS  q  at  5.60  bar  for  R-13b1  (Engelhorn) 


Logarithmic  Curve  Fit 

Pts  Plotted  =  7  Offscale  Pts  =  0 

Regression  Equation: 

Y  =-20.435+  (3.15063  )  LNX 

Correlation  Coefficient  =  .984194 
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Pressure  vs  Departure  Diameter  for  R13B1 
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